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Abstract 

The aim of this paper is to show that one can achieve convergence rates of N~ a+S 
for a > 1/2 (and for 5 > arbitrarily small) for nonperiodic a-smooth cosine series 
using lattice rules without random shifting. The smoothness of the functions can 
be measured by the decay rate of the cosine coefficients. For a specific choice of the 
parameters the cosine series space coincides with the unanchored Sobolev space of 
^ smoothness 1. 

— We study the embeddings of various reproducing kernel Hilbert spaces and nu- 

merical integration in the cosine series function space and show that by applying the 
so-called tent transformation to a lattice rule one can achieve the (almost) optimal 
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rate of convergence of the integration error. The same holds true for symmetrized 
lattice rules for the tensor product of the direct sum of the Korobov space and 
cosine series space, but with a stronger dependence on the dimension in this case. 

1 Introduction 

Quasi-Monte Carlo (QMC) rules are equal weight quadrature rules 

n=0 



which can be used to approximate integrals of the form 

/ f( x ) dx i 

see [10, 23, 28] for more information. In QMC rules, the quadrature points {x , xi, . . . , x N _i} 
are chosen according to some deterministic algorithm. One can show that the convergence 
rate of the integration error J2n=o f( x n) — Jj ]> f( x ) da? depends on the smoothness 
of the integrand and some property of the quadrature points. 

There are several deterministic construction methods for the quadrature points. One 
such method yields so-called digital nets. These yield a convergence of the integration 
error of iV -1 (log iV) s for functions of bounded variation [10, 23]. Higher order digital nets, 
using an interlacing factor of d, yield a convergence rate of order jV~ min ( a ^ (log N) sa for 
integrands with square integrable partial mixed derivatives of order a in each variable 
[8, 9]. 
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An alternative to digital nets are lattice rules [23, 28]. In this case, for a given positive 
integer A" one chooses a vector g £ {1, . . . , N — 1} S called the generating vector and 
defines the quadrature points by 



Here, for a real number x, the braces {x} denote the fractional part, i.e., modulo 1. For 
vectors, the fractional part is taken component wise. It is well known that there are 
generating vectors for lattice rules for which the integration error converges with order 
jy-a+<5 ^ > q^j f Qr sm oothness a > 1, but with the restriction that / and its partial 
derivatives up to order a — 1 in each variable have to be periodic [19, Theorem 18, p. 120] 
and [29, 30]. Fast computer search algorithms for such vectors are known from [25, 26]. 
Hence, in order to be able to benefit from the fast rate of convergence, one needs to apply a 
transformation which makes the integrand (and its partial derivatives) periodic. This can 
cause some problems though and is not always recommended [19, 22]. Since arbitrarily 
high rates of convergence can be obtained using digital nets for nonperiodic functions, 
the question arises whether this is also possible for lattice rules. Until now, lattice rules 
achieve a convergence rate of at most N~ 1+s (5 > 0) for nonperiodic integrands (via 
estimates of the star-discrepancy [23]). If one applies the so-called tent transformation 
and a random shift, this rate of convergence can be improved to N~ 2+s for any 5 > for 
the worst-case error in an unanchored Sobolev space of smoothness 2, see [15]. 

In this paper we present quadrature rules which achieve a convergence rate of order 
N~ a+S , 5 > 0, for nonperiodic functions with smoothness a > 1/2. The way we measure 
smoothness in this paper is slightly different from the setting used in, for instance, higher 
order digital nets [8, 9] though. We consider functions / which can be represented by a 
cosine series. Note that every continuous function / £ L 2 ([0, 1]) can be represented by a 
cosine series (see [18, Theorem 1] for this basis over [—1,1]). This follows from the fact 
that the functions 



are L 2 -orthogonal and complete. 

As mentioned above, the "smoothness" of the cosine series in our context is measured 
by the rate of decay of the cosine coefficients. To illustrate this, consider a one-dimensional 
function / : [0, 1] — > R given by its cosine series 




for < n < N. 



cos(7rx), \pl cos(7r2x), \/2 cos(7r3x), . . . 



oo 



/co S (0) + ^2 f cos (k)V2cos(-Kkx) 



k=l 



oo 



oo 



/(0) + fco S (2k)V2cos(n2kx) + 



cos 



(2k - ljv^cos^ k - l)a?) 



fe=l k=l 
C ~\~ fpeiipC) ,/nonper (*^J j 



where 




for k — 0, 



for k £ N. 



(1) 
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The sum over the even frequencies is a 1-periodic function f per over [0, 1]. If the coeffi- 
cients f(h) decay with order h~ a then / per is a-times differentiable in the classical sense. 
However, this does not apply to the sum over the odd coefficients. For instance, the cosine 



and hence the odd coefficients converge with order k~ 2 only, although x is infinitely times 
differentiable. 

Below we introduce a reproducing kernel Hilbert space based on cosine series, with the 
smoothness measured by the decay rate of the cosine coefficients. Although the smooth- 
ness of a cosine series measured by the differentiability of the series can be larger than the 
decay rate of the cosine coefficients suggests, the opposite can not happen. That is, we 
show that the reproducing kernel Hilbert space based on cosine series is embedded in the 
unanchored Sobolev space with the same value of the smoothness parameter. The case 
of smoothness 1 provides an exception, since there the cosine series space and the unan- 
chored Sobolev space coincide. Various reproducing kernel Hilbert spaces are introduced 
in Section 2 and their embeddings are studied in Section 3. 

In this paper we present two methods which allow us to achieve a higher convergence 
rate for smoother nonperiodic functions using lattice rules, namely: 

1. application of the tent transformation to the integration nodes; 

2. symmetrization of the integration nodes. 
The tent transformation, <p '■ [0, 1] — > [0, 1], 



is a Lebesgue measure preserving function. The idea of using this transformation in con- 
junction with a random shift for integration based on lattice rules comes from Hickernell 
[15] and was also used in [6] for digital nets. In contrast to these works, here we do not 
rely on a random element in our quadrature rules. We show that a tent transformed 
lattice rule achieves an integration error of order jV -Q+<s , for any 5 > 0, for functions 
belonging to a certain reproducing kernel Hilbert space of cosine series with smoothness 
parameter a. This result follows by showing that the worst-case error in the cosine series 
space for a tent-transformed lattice rule is the same as the worst-case error in a Korobov 
space of smooth periodic functions using lattice rules. Thus all the results for integration 
in Korobov spaces using lattice rules [7, 12, 20, 25, 26] also apply for integration in the 
cosine series space using tent-transformed lattice rules. In particular for smoothness 1, 
this yields deterministic point sets for numerical integration in unanchored Sobolev spaces 
with the same tractability properties as for numerical integration in the Korobov space. 

Furthermore, we also use symmetrized lattice rules. We show that these rules achieve 
the optimal order of convergence for integration of sums of products of cosine series and 
Fourier series. We apply the transformation x t— > 1 — x to each possible set of coordinates 
separately, so that if we start off with N points we get 0(2 S ~ 1 N) points (see Section 4.2). 
This symmetrization approach is also mentioned in [19, 28, 34] and is one of the symmetry 



series for £>—)•£ — - 



is given by 




(2) 
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<j>{x) := 1 - \2x - 1 
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groups applied in the construction of cubature formulae, see, e.g., [4, 13]. We prove that 
a lattice rule symmetrized this way achieves an integration error of order N~ a+S , 5 > 0, 
for functions belonging to a certain reproducing kernel Hilbert space of cosine series and 
Fourier series with smoothness parameter a. The advantage of symmetrized lattice rules is 
that functions of the form cos(7r(2A; — 1)x), where k is a nonnegative integer, are integrated 
exactly and hence only the smoothness of the periodic part determines the convergence 
rate. Note that the decay rate of the cosine series coefficients of the periodic part and 
of the Fourier series coefficients part coincides with the classical smoothness. Thus the 
problem with functions where the smoothness in terms of differentiability differs from the 
rate of decay of the cosine series is overcome using symmetrization. For instance, the 
function x h-> x is integrated exactly using symmetrization. However, a disadvantage of 
the symmetrization compared to the tent transformation is that the number of function 
evaluations grows exponential in the dimension and therefore symmetrization is only useful 
in smaller dimensions. 

For both methods, the rates of convergence we obtain are essentially optimal by an 
adaption of the lower bound of Bakhvalov [3], which is presented in Section 4.3. 

In the next section we introduce four reproducing kernel Hilbert spaces, the unan- 
chored Sobolev space, the Korobov space, the cosine series space and the sum of the 
cosine and Korobov space. Since the unanchored Sobolev space and the Korobov space 
are frequently studied in the literature, we study the relations among these four spaces 
in Section 3 to put our results into context. It is shown that the Korobov space and 
the cosine series space differ, but both are embedded in the sum of the cosine series and 
Korobov space, which is itself embedded in the unanchored Sobolev space. In Section 4 
we study numerical integration in the cosine series space using tent-transformed lattice 
rules and numerical integration in the sum of the cosine series and Korobov space using 
symmetrized lattice rules. Numerical results are presented in Section 5 and a conclusion 
is presented in Section 6. 

We write Z for the set of integers, N := {1,2,...} for the set of positive integers and 
No := {0, 1, 2, . . .} for the set of nonnegative integers. We also write W + := {x G R : x > 
0} s . Furthermore, for s G N we write [s] := {1, . . . , s}. 

2 Reproducing kernel Hilbert spaces 

In this section we introduce several reproducing kernel Hilbert spaces [2] . For a reproduc- 
ing kernel K : [0, 1] x [0, 1] — > R we denote by TL(K) the corresponding reproducing kernel 
Hilbert space with inner product (•, -)k and corresponding norm ||/||k = a/ (/, f)x- For 
any y G [0, 1] we have K(-,y) G H{K) and we have the reproducing property 

f(y) = (f,K(-,y)) for all y G [0,1] and / G H{K). 

Further, the function K is symmetric in its arguments and positive semi-definite. 

For higher dimensions s > 1 we consider tensor product spaces. The reproducing 
kernel is in this case given by 

s 

K s {x,y) = Y[K(xj,yj), 
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where ) and y = (yi, . . . ,y s ). Again, the corresponding reproducing kernel 

Hilbert space is denoted by T-L(K S ), the corresponding inner product is denoted by (•, -)k s 
and the corresponding norm by ||/||k- s = a/ (/, J)k s - 

For further information on reproducing kernel Hilbert spaces we refer to [2] (or [10, 
Chapter 2] for reproducing kernel Hilbert spaces in the context of numerical integration). 

2.1 The unanchored Sobolev space 

The unanchored Sobolev space is a reproducing kernel Hilbert space l-L(Kf^) with repro- 
ducing kernel Kf^ : [0, 1] x [0, 1] — > R given by 

K^(x,y) := l+jB 1 (x)B l (y)+j B ^ X ~ yl \ 

where Bi(z) = z — 1/2 and B 2 (z) = z 2 — z + 1/6 are Bernoulli polynomials and 7 > is 
a real number. The inner product in this space is given by 

1 rl l rl 



(f,g)Ki° b = f{x)dx g(x)dx + - f(x)g'(x)dx. 
' 7 Jo Jo 7 Jo 

In higher dimensions s > 1 we consider the tensor product space "H(i^^)®- • •®'H{Ky? s ). 
The reproducing kernel is in this case given by 

s 

3=1 

where x = (xi, . . . ,x s ), y = (yi, . . . , y s ) and 7 = . . . , 7s ) g R s + . 

The reproducing kernel Hilbert space T-L^Kf^) can be generalized to higher order 
smoothness. For a G N, consider the reproducing kernel : [0, 1] x [0, 1] — > R given 

by 

K^ix y)-=l + iY B ^ B ^V) _ { _ l)a B 2a {\x-y\) 

where B T is the Bernoulli polynomial of order r. The inner product for this space is given 
by 

(f,g)K° a °*= f(x)dx g(x)dx 
'0 Jo 



+ -Y^ I f [T) (x)dx I g ir) (x)dx + - I f a \x)g {a \x)dx. 
7^io Jo 1 Jo 



To obtain reproducing kernel Hilbert spaces for the domain [0, l] s , we consider again the 
tensor product of the one-dimensional spaces. This space has the reproducing kernel 

s 

3=1 

Quasi-Monte Carlo rules in H(K^ s ) which yield the optimal rate of convergence of 
order N~ a+S for any 5 > where studied in [9]. 
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2.2 The Korobov space 

For a > 1/2, ft£Z and 7 > we define 



r*a{h) { 7 | /i |_ 2Q if h _^ ] (3) 



For fa, = (hi, . . . , h s ) eZ s and 7 = (71, . . . ,7 S ) € Ei, we set 



:= JJ'W^)- 



The Korobov space is a reproducing kernel Hilbert space of Fourier series. The repro- 
ducing kernel for this space is given by 



and in higher dimensions s > 1 by 



y) 



j=i hez s 

where the "•" denotes the usual inner product in W. Let the Fourier coefficient for a 
function / : [0, l] s — > R be given by 



f(h) := f f(x)e^ h - x dx. 



'[0,l] s 

Then the inner product in the reproducing kernel Hilbert space H(K^ s ) is given by 

(f,9)K* %s =J2f( h ^ r «n,s( h )- 



h&L s 



The corresponding norm is defined by H/Hxta, = J(f,f)gg* y 
2.3 The half-period cosine space 

The half-period cosine space is a reproducing kernel Hilbert space of (half-period) cosine 
series with reproducing kernel 

00 

K™(x,y) := 1 + r aa (k) V2 cos(irkx) a/2 cos(7r%) , 

k=l 

where r Qj7 is defined as in (3) and where a > 1/2 and 7 > 0. The inner product is given 
by 

00 

k=0 
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where / cos and g cos are the cosine coefficients of / and g, respectively, as defined in (1). 

We can generalize the reproducing kernel Hilbert space *H(K%t) to the domain [0, l] s 
by setting 

s 

K:;:;jx.y): ]jK;;::(.r r!Jj ). 

where 7 = (71, ...,7,)£K' and x = {x u ...,x s ),y = (yi,...,y s ) G [0, 1} S . The inner 
product is then given by 

(/>9)^, s = Yl /cos(fc)^cos(fc)r^ 1 7 S (fc), 
feeNg 

where the multi-dimensional cosine coefficients for a function / : [0, l] s — > R are given by 

Jcos(fc) := / f(x) 2l fc l°/ 2 TT cos(7r%^) daj, 
where for k = (ki, . . . , k s ) G Nq we define |fc|o := \{j G [s] : kj 7^ 0}| to be the number of 



nonzero components in k. The corresponding norm is defined by ||/||jf«» s = \/ (/, f)x\ 
in particular we have 

|7cos(fc)| 2 _ ^ o _ W0 l7cos(l^l)| 2 



COH 



where \h\ = (\hx\, . . . , \h d \). 



2.4 The sum of the Korobov space and the half-period cosine 
space 

In this section we introduce the kernel 

1 1 00 

= 2 2 r ^ h) e 2niHx ~ y) + g + E r «^) cos(vrA;a;) cos(vr%) 

hel k=l 



1 + 7 k 2a (cos(2iik(x — y)) + cos(7rA;x) cos(nky)) . 



k=l 



The space resulting from the sum of kernels is studied in [2, Part I, Section 6]. The norm 
in the reproducing kernel Hilbert space 'H(K^ +COB ) is then defined by 



K; 



mrn ^ ( llAorll Jfkor + 1 1 /cos 1 1 lc co; 



kor+cos 

a 'T /=/kor + /c, 

kor"i 



where the minimum is taken over all functions /k or G 7i(K%£) and / cos G %{K^) such 

that / = /kor + /cos- 
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For dimensions s > 1 we define the reproducing kernel by 

K%*T(x, y) := n \ Vi) + K% fo, %)) • 

Thus the space n(K^+ cos ) is the tensor product U{K^+ cos ) <g> ■ ■ ■ <g> ^(^° 7 r + cos ). For 
u C [s] we define 

^ssn*. ?/) ; = 2- s n «•'>//,) n 

ieu je[s]\u 

where as usual an empty product is considered to be one. This is a reproducing kernel 
for the space 



' 1 7>^COS * 

jeu J V'e[s]\u 



with the inner product 

(/) 9)ku°^ + sZ S = 2 S /u,kor+cos(^U) fc[s]\u) 5 , u,kor+cos(^ i U) fc[s]\u) ^^sl^U) fc[s]\u)i 

^u6Zlu| fcwXueN J-|u| 

where 



/u,kor+cos(^u,fc [s] \ u ) := / f(x)(V2f^Tle- 2 ^ 11 cos^A^d*. 

J[o,i]s tt 

Clearly we have that 



jeu je[s]\u 



uC[s] 

3 Embeddings 

In this section we investigate the relationships between the spaces introduced above. 

For two reproducing kernel Hilbert spaces 'H(Ki) and %{K 2 ) we say that 'H(Ki) is 
continuously embedded in T-L(K 2 ) if 

and if 



k<C||/|k for all / G %{Ki) 
for some constant C > independent of /. We write 

in this case. 

On the other hand it is possible that T-L(Ki) is not a subset of H(K 2 ), i.e., there is a 
function in l-L{Ki) which is not in %{K 2 ). In this case we write T-L{K]) ^ H(K 2 ). 
If ft(tf 2 ) and n(K 2 ) ^ H^) we write 

U{K X )^U{K 2 ). 
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3.1 The Korobov space and the unanchored Sobolev space 

It is well known, see, e.g., [24, Appendix A], that the Korobov space is continuously 
embedded in the unanchored Sobolev space 

Conversely, as is also well known, for instance the function g : [0, l] s — > IR, g{x) = x\ is 
in %{K S ^ S ) for all a£N, but not in 1-L(K^*), since g is not periodic. Thus 

We note that for / 6 H(K^ S ), a 6 N, we have that ||/||#-sob = ||/||x sob where 

' a '^ ,s a,'y(27r) — ^ a ,s 

j(2ii)~ 2a denotes the rescaled sequence (7i(27t)~ 2q , . . . , 7 s (27r) _2a ). 

3.2 The half-period cosine space and the unanchored Sobolev 
space 

We now consider the half-period cosine space and the unanchored Sobolev space. For 
a = 1 there is a peculiarity. 

Lemma 1. We have 

K$(x,y) = 1 + -1 ]T -2cos(fcvrx) cos^vn/) = K^. 2 (x,y). 

k=l 

For completeness we include a short proof. 



Proof. In the following we calculate the cosine coefficients of Kf 3 }]. It is easy to check 
that we have J* Kf°^(x, y) cos(nnx) cos(iTmy) dx dy = for (n,m) 6 {(0,/c) : k > 
0} U {(k, 0) : k> 0}. Further we have K$(x, y) dx dy = 1. 



We have 



1 .__ p 2->rik(x-y) 

B 2 {\x-y\)=B 2 {{x-y}) = — 2 £ 

fcez\{o} 



E X1 cos(2nkx) cos(2irky) sin(27rfca;) s'm(27rky) 

7T 2 k 2 ^ ix 2 k 2 

k=l k=l 



Using (2) we obtain 



(* -i)(v-i) + B2 K vl) = £ n * i2k _ i) 6 2(2/ _ 1)2 co « 2k ~ ^ cos ^ 21 - x )y) 



k 

oo 



Ecos(27rA;x) cos(27r%) ^-^ sin(27rA;x) sin(27r/o/) 
2v 2 k 2 + ^ 2n 2 k 2 

k=l k=l 



This immediately implies that 



1 r l 

Kf^V2 cos(irmx) cos{imy) dxdy = 



o jo 



if m is even and n is odd, or m is odd and n is even, or m, n are even with m ^ n. If 
m = n = k for even k > 0, we obtain 

Kl°^(x, y)V2 cos(irkx)\/2 cos(nky) dx dy 



o jo 



7T 2 k 2 

Now let m, n > be odd. We have sm(2nkx) cos(irmx) dx = ^fe^z^S) and 



therefore 

f 1 f 1 8 16 

/ / Kl°hx,y)V2cos(irmx)V2cos(Trny)dxdy = - + — 



' (4k 2 - m 2 ) (4k 2 - n 2 )' 



For m ^ n we have (4fc2 _ m2) 1 (4fc2 _ n2) = and further ££° =1 psj^p = ^f^. 

Thus we obtain 

/ / Kl^(x,y)V2cos(Tcmx)V2cos(Tcny)dxdy = < 7 . f = ' 

Note that the cosine series for Kl°z converges absolutely. Since the function Kf 3 ^ is 
continuous, the cosine series converges to the function pointwise. This completes the 
proof. □ 

The above lemma and Mercer's theorem also yield the eigenfunctions of the operator 

T(g)(y)= f K$(x,y)g(x)dx. 



These are 1, a/2 cos(7rx), a/2 cos(27ra;), \/2cos(3tcx), . . . and the corresponding eigenvalues 
are l,7T- 2 ,(7r2)- 2 ,(7r3)- 2 ,.... 

Remark 1. An analoguous result for a slightly different reproducing kernel Hilbert space 
was established in [33]. For this space one obtains the same set of eigenfunctions. For 
the anchored Sobolev space the eigenfunctions are slightly different and have been found 
in [32]. 

We thus find 

U(K^) = U(K™^) and ll/lljq* Hl/ll*^_ a for all / G H(K^). 
The same also applies for the higher dimensional tensor product space 

H(Ki°* s ) = U{K™^ S ) and \\f\\ K ^ = \\f\\ K ^_ 2s for all / e H(K^ S ), 

where jn~ 2 denotes the sequence (7i7r~ 2 , . . . ,7 s 7r~ 2 ). Thus 

Now consider a > 1. The function x i— >■ x belongs to H(K S ^) for all a G N. On the 
other hand we have 

X= 2~V 2 ^ k 2 ' 

k=l 
k odd 
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Hence the function x (->■ x is not in H(K^) for a > 3/2 and therefore 

n(K^)tn(K^) foraeN,a>2. 
Conversely, let a G N, a > 2. Let / G H(K™*) be given by 

oo 

f(x) = f cos (0) + ^^ os (A;)v / 2cos(7rA;a;) 
fc=i 

with 



|7cos(0)| 2 + -V|7 C os(A:)| 2 |A;| 2a <oo. 
Then for 1 < r < o we have 

oo 

f {T \x) = Y,fcosm-iy T/2 XknyV2<j )T (*kx), 



k=l 



where <f> T (z) = cos(z) for r even and <fi T {z) = sin(^) for r odd and where for a real number 
x, \x] denotes the smallest integer bigger or equal to x. Thus 



1 

7 



/ f (T) (x)dx <-/ \f^\ X )fdx=—J2\fcos(k)^<^^ 

Jo 7 Jo 7 ^ 



/("COS • 

Q >7 



Thus 



1/2 / 7r 2(a+l)_ 1 N V2 



I> 2r 

,r=0 

Thus we have 

This result can be generalized to the tensor product space, thus 

3.3 The half-period cosine space and the Korobov space 

For a = 1 the embedding results for the half-period cosine space and the Korobov space 
follow from the previous two subsections. 
Let now a > 1. Let 

f(x) = sin(27r:r). 

Then / G T-L(K^) for all a > 1/2. On the other hand we have for k G N 



fcos{k) 

Thus 



2 



/ sin(27rx)v^cos(7r/i;x)dx = < ^( 4 - fc2 ) lf k 1S 
Jo otherwise. 



1 V- ,7 /m,2,2«_ !32^ (2fc - 1) 2 « 



*™ ~7 S l/c ° s(A;)|2A;2a ~7 7r2 U((2^-l) 2 -4) 2 



fe=i 

fe odd 



11 



Thus we have 



oo for a > 3/2. Thus 




hor{h)= [ cos(7ra:)e- 2 ^da ; = 
Jo 




For any a > 1/2 we have 



11/11 



1 16 



2a: 



h£Z\{0} 



(4/z 2 - 1) 



which implies that / ^ T~L{K ( 



"kor 

a, 7 



As we need a > 1/2 we therefore have 



3.4 Embeddings of the sum of the Korobov and half-period co- 
sine space 

Since K*°* +cos = ±K%% + \K™ we obtain from results from [2] that 



On the other hand, as it was shown above, we have H(K%*), 7i(K^t) c — >■ 'H(K^), thus 

We have seen above that for a > 1/2 we have H(K™*) (£ H(K^), thus 

and for a > 3/2 we have U(K^) <£ U{K C °*) and therefore 

U{K^ +cos ) £ U{K™) for a > 3/2. 
For a = 1 we have T-L{Kf^) <=> 'H(K^) and therefore 



"kor+cos 
a, 7 

"kor+cos 
ct,7 
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3.5 Summary of embeddings 

We summarize the obtained embedding results in the following theorem. 
Theorem 1. For a = 1 we have 

and 

For a G N with a > 1 we have 
and for a G K we have 

H(K%*r) t H(K^ S ) for a > 3/2, 
U(K^ S )<£H(K™ >S ) for a > 3/2. 

We do not know whether U{K s a °^ s ) differs from n(K^+ COB ) for a > 1. 



4 Numerical integration 

We now study the worst-case error for QMC integration. As quality measure for the QMC 
algorithm we use the worst-case integration error. 

Let P = {xq, . . . , a?Ar_i} and let H{K) be an arbitrary reproducing kernel Hilbert 
space with reproducing kernel K and norm || ■ Then the worst-case error for QMC 
integration in 'H(K) using the point set P is defined as 



e(H(K); P) := sup 
fen(K) 

ll/lk<l 



1 

/ f(x) dx - — V /(av 



see, e.g., [10, 24] for a general reference. We use the following formula for the square 
worst-case error (see [10, Proposition 2.11] or [14]): 



N-1 „ , N-1 



e\U{K)- P) = I K(x, y) dx dy - ^ V / K(x, x n ) dx + -L V a;„/). 

Integration in the Sobolev space 'H{K S ^ ( S ) has been considered in [9, 10] and integra- 
tion in the Korobov space has been studied for instance in [7, 8, 17, 19, 20, 25, 26], as well 
as other papers. In this paper we study numerical integration in the half-period cosine 
space and in the sum of the Korobov space and the half-period cosine space. For the for- 
mer space we use tent-transformed lattice rules and for the latter one we use symmetrized 
lattice rules. 
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4.1 Numerical integration in the half-period cosine space 

We now study numerical integration in the half-period cosine space using tent-transformed 
lattice rules. For a nonnegative real number x we denote the fractional part of x by 
{x} = x — [x\. For a vector x of nonnegative real numbers, the expression {x} denotes 
the vector of fractional parts. A lattice point set with N > 2 points and generating vector 
g G {1, . . . , iV — 1} S is given by 

P(9,N):={{^r}:0<n<N}. (4) 

For x G [0, 1] we define the tent-transformation by <f>(x) — 1 — \2x — 1| and for vectors 
we apply the function <fi component-wise. The tent-transformed lattice point set is now 
given by 

PMN):= {<t>({^r }) :0<n<N}. 

We call a lattice rule which is based on P^,(g, N) a tent-transformed lattice rule. 

The following theorem gives a useful formula for the worst-case integration error in 
HiK^ s) of tent-transformed lattice rules. 

Theorem 2. The squared worst-case error for QMC integration in the half-period cosine 
space H(K™^ iS ) using a tent-transformed lattice rule is given by 

heL±\{o} 

where L 1 - := {h G Z s : h ■ g = (mod N)} is the dual lattice. 
Proof. Let / G 7~L(K^ S ) with ||/||^cos s < oo and with expansion 

s 

fix) = 7cos(fc)(v / 2) |fc| ° Hcosinhxi). (5) 
ke% j=i 

For any k G N we have 

cos(7r k(f){x)) = cos{2itkx) for all x G [0, 1], 

and hence 

k&% j=i 
= ^(V^)"*l«»/„(fc)n 008(2^-^) 

ke% j=i 
= J2(^y lhl "f^(\h\)e 2mn ^ 9)/N . 

Therefore we obtain 

N-i f / , N-l \ 

h £' {* f ) - [01] , m ^ - E fc/W-<w) £ E ^ )IN ) ■ 

n=0 " / l U ' 1 J 0^h.eZ s \ n=0 / 
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The sum in the braces is a character sum over the group Z/iVZ which is one if h ■ g is a 
multiple of N and zero otherwise. From this we get 



n=0 IU h£L±\{0} 



(6) 



From this formula and an application of the Cauchy-Schwarz inequality we obtain 

JV-l 

N 



Y r ^A h ) 

heL x \{0} 



1/2 (V2)-^ f cos (\h\) 



1/2 



<( E r m {h)\ [Y 

1/2 



2-' fa '"|/cos(|fa| 



1/2 



E 

v fc£lA\{0} 



I/cos 
a,-)',, 



Here we obtain equality by (6) for the function with cosine series coefficients given by 

fcos(k) 



(V2) |fe|o r Q)7)S (fc) for k e {\h\ : h e L 1 \ {0}} 
otherwise. 



Hence the result follows by the definition of the worst-case error. 



□ 



The above result shows that the square worst-case error of tent-transformed lattice 
rules for numerical integration in the cosine space coincides with the square worst-case 
error for numerical integration in a Korobov space using the same lattice rules but without 
applying the tent-transformation, that is, 

2 



e 2 (n(K^ >s );P(g,N)) = -1 + £ r an , s (h) 



1 N ~ l 
N ^ 



n=0 



E 



(h) 



heL-L\{0} 

= e 2 (H(K™ tS )-P^g,N)). 

Thus all the results which hold for the worst-case error in the Korobov space using lat- 
tice rules, also hold for the worst-case error in the half-period cosine space using tent- 
transformed lattice rules. This applies for instance to the component-by-component con- 
struction [7, 20, 30, 29] and fast component-by-component construction [25, 26], general 
weights [12, 21] and extensible lattice rules [5, 11, 16, 17]. 

Corollary 1. Using the fast component-by- component algorithm one can obtain a gener- 
ating vector g e {1, . . . , iV — 1} S such that 



e{U{K™ >s )- P^g, iV)) < C Q , 7 , s , T (iV - 1) 



-t/2 
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for all 1 < t < 2a, where the constant C a ~ ST > is given by 



-/ 2 , x r/2 



= ( £ 7y r (2C(2«/r))l"l = -1 + H(l + 2C(2«/r) 7 ] / H 

0^uC[s] / V 3=1 



Note that for certain choices of weights 7, the upper bound can be made independent 
of the dimension s and we then obtain (strong) tractability results. See [7, 12, 20, 24] for 
a discussion of tractability results which apply in this context. 

Since the cosine series space H(Ki°~ s ) and the unanchored Sobolev space 7~L(K^ S ) 
coincide, we also get the following result. 

Corollary 2. Using the fast component-by- component algorithm one can obtain a gener- 
ating vector g G {1, . . . , N — 1} S such that 

e(n(K^ >s ),P^g,N)) < C ljW (iV - l)^ 2 , 

for all 1 < t < 2, where 

Ci„,. lT = ( £ 7u 1/T (2C(2/r))l u ^ = (-1 + ;Q(1 + 2C(2/r) 7 f) X 

\0^uC[ s ] / V 3=1 

Remark 2. An analoguous result can also be obtained for the space considered in [33], 
since the eigenfunctions are the same as for H{Kf^). 

This result is a deterministic version of the main results in [7, 20], where a random 
shift was required to achieve this bound. The tractability results of [7, 12, 20] also apply. 

4.2 Numerical integration in the Korobov plus half-period co- 
sine space 

We now study numerical integration in the space 'H(K^^ s cos ) using symmetrized lattice 
rules. Let ) and let u C [s]. Then let sym u (:r) denote the vector whose jth 

coordinate is Xj if j G u and 1 — Xj otherwise, i.e., sym u (x) = (y\, . . . ,y s ) with 



y.i 



1 — Xj if j G u, 
Xj if j ^ u. 



For a lattice point set as in (4) let 

Ps ym (<7, AO := {sym u ({^ }) : < n < N,u C [s]} . 

We call a lattice rule which is based on P S ym(fl r ) N) a symmetrized lattice rule. Note that 
P sym (g, N) consists of 0(2 s ~ l N) elements as we show next. 

Lemma 2. The number of nodes in the symmetrized lattice rule P sy m{9,N) is given by 
2 S - 1 (N + 1) if 2 \ N and 2 S ~ 1 N + 1 if 2 | N. 
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Proof. The argument comes from [34]. We have the following symmetry 

kgj =N-(N- k)g 3 (mod N), for all < k < N and j = 1, . . . , s, 

which corresponds exactly to x^j = 1 — x^-k ,j- This means we only have to evaluate and 
symmetrize the points for < k < N/2. 

1. For < k < N/2 we have 2 s (iV - 1) /2 points if 2 f N and 2 s (iV/2 - 1) if 2 | N. 

2. For k = symmetrization returns all 2 s corner points. 

3. If 2 | N then for k = N/2 we have Xn/2 = (§, • • • , ^) and thus no symmetrization is 
needed. 

Counting the number of function evaluations now gives the result from above. □ 

The analysis used in the previous proof can also be used in an implementation where 
further stream lining can be done by noticing that x H- 1 — x is its own inverse and thus 
the 2 s symmetric points for < k < N/2 can be constructed easily in gray code ordering. 
Nevertheless, we increase the number of points by a factor of 2 S_1 and so this is only 
feasible for moderate dimensions. For the derivations below we symmetrize all N points 
for ease of notation. 

The following theorem gives a useful formula for the worst-case integration error in 
%(ifk°r+ cos ) f symmetrized lattice rules. 

Theorem 3. The squared worst-case error for QMC integration in the sum of the half- 
period cosine space and the Korobov space 'H(K^^ s cos ) using a symmetrized lattice rule is 
given by 

e 2 (n(K^ry,PsyUg,N))= £ r^h), 

/ie£ x \{o} 

where L 1 - := {h G V : h ■ g = (mod iV)} is the dual lattice. 

Proof. Let / G ^(^+ cos ) with \\f\\ K ^r+cos < oo. Let C [s], then by [2, Part I, Sec- 
tion 6] there are functions / G H(K%;+™) = (<g>. £0 U{\K^)) ® (® i6MV ,«(§ff, 
given by 



a 



1i- 



h ezi»i fe wv ,eN*- 101 J 6D i6[«]\B 



kor+cos 7 1 1 J X) 1 1 iv'kor+cos • 



such that 

Eli/- 

tC[s 

Note that for fceZ and x G M, 



cos(7r/ca;) + cos(7r/c(l — x)) 



2 cos(7rfcx) = e 71 " 1 * 131 + e mkx if A; is even, 
if k is odd. 
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For given t) C [s] we therefore have 

N-l 

1 

2~ 



1 JV_1 n 

^EE/»h({f})) 



n=0 uC[s] 

AT-1 



1 JV_± ~ 

n=o ht>eZ M fcwXoeNr i»i 

J] (e 2nih ^ /N + e -^ jngj /N^ "Q ( cos ( 7r A; J {n^/iV}) + cos(7r%(l - {n^/iV}))) 



x 

je» je[s]\v 
1 iV - 1 

^EE E (v^) |fe '^" l0 / B ,ko r+ cos(^,2fc [sAt ,) 

n=0 ht>eZ M fc[s]Xo6Nr H 
x (^Kihjngj/N _|_ ^-2-Kihjngj /N ^ ^wikjugj /N _|_ ^-2-Kikjngj /N^ 

jet> j6[s]\o 

iV-1 

2 S 7V 



^ ^ ^ 2l B l(v^) |h w\» lo 2 s -H-^wv»IV B ,ko r+ cos(^,2|h [s] \ D |)e 2 -^^ 

n=0 freZ 3 



he1 s \ n=0 



The sum in the braces is a character sum over the group Z/iVZ which is one if h ■ g is a 
multiple of TV and zero otherwise. From this we get 



rj,=n „rr*i "M - 1 ! 



n=0 uC[s] 



^ (v^)-l h w\«-l%kar+oo B (Ai B ,2|h [a] \ B |). (7) 
hei/ J -\{o} 



Thus we find 
1 

2~ 



5vEE/hO))-/ pi /(^ 

n=O u C[s] ^l ' 1 ! 
^ E E (^"'^'^Akor+coe^O^lh^l)! 

bC[s] heL x \{o} 

V- V- r ,u/ 2 (v^)~ |,t '^'' l0 |7 ,kor+co S (fcD,2|/t[ fl] \ t ,| 

bC[ s ] heL-L\{o} a ' 7 ' sv 7 

/ \ 1/2 / ~ \ 1/2 

< (E E E E r "" , " l/ -ft' %J " 

/ \ 1/2 / _ ^ 12 

< 2' £ r TO (fc) !^E 2S E E IWos( ^' fc[ ^ )|2 



heL x \{o} 
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1/2 / \ 1/2 

2 



^-kor+cos 



< F E 

^£L-L\{0} / \ BC[ S ] 

\ 1/2 

^ ( E ^.-CO 
^£L-L\{0} 

Thus the result follows. □ 

Again we can relate the worst-case error e 2 {T-L(K^ + s cos ); P sym (g, N)) to the worst-case 
error in a Korobov space. We have 

e(U(Kl™y° s ); P sym (g, N)) = e(U(K^ s ); P(g, N)). 

Thus the results for integration in the sum of the half-period cosine space and the Korobov 
space using symmetrized lattice rules are the same as in a Korobov space using lattice 
rules. In particular, the component-by-component algorithm can be used [7, 20, 30, 29] 
and also its fast version [25, 26], general weights [12, 21] and extensible lattice rules 
[5, 11, 16, 17]. 

Corollary 3. Using the fast component-by- component algorithm one can obtain a gener- 
ating vector g G {1, . . . , iV — 1} S such that 

e(H(K^ry,P S yU9,N)) < C a ^ T (N-l)-^ } 

for all 1 < r < 2a, where the constant C antSyT > is given by 

Co, w = ( E lt/ T (n^/r))A = (-1 + f[(l + 2C(2«/r) 7 ] /T ) X 

\0^uC[s] / V j=l 

Due to the fact that the number of points of P sym (g, N) is M = 0(2 S ^ 1 N) we do not 
get tractability results. Notice that in terms of the number of points one gets (N— l)~ r / 2 « 
2( S ~ 1 ') T / 2 M~ T / 2 , which also implies a strong dependence on the dimension. 

A consequence of the symmetrization procedure is that all functions of the form 

s 

E b k U :.,k a JJcos^/cjX,) for all b ku „. )ks G R 

fei,...,fc s eN j=i 

ki,...,k s odd 

are integrated exactly. Likewise, all polynomials of the form 

s 

E a ku-,ks ~ l / 2 P f o r all a feli ... )fca G R 

fci,...,fc s eN j=i 

ki,...,k s odd 

are integrated exactly. This is because all the odd frequencies in a cosine series are 
integrated exactly. Specifically for the half-period cosine space we can state the following 
result. 
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Corollary 4. The squared worst-case error for QMC integration in the half-period cosine 
space H(K™ s ) using a symmetrized lattice rule is given by 

e 2 (U(K™ ;S );P sym (g,N)) = £ r w (2fc), 

heL x \{o} 



where L 1 - := {h G V : ft. ■ g = (mod iV)} zs £/ie ditaZ lattice. 
Proof. Similar to the proof of Theorem 3 with = 0. 



□ 



4.3 A lower bound on the worst-case error 

We prove the following lower bound for integration in the half-period cosine space. Let 
P = {x , . . . , a?AT-i} C [0, l] s be an N element point set and let w = (u> , • • • , wjv-i) be 
an arbitrary real tuple. Let 



e(H(K; 



sup 

feH(K°% s ) 



» N-l 

/ f{x) dx - V w n f{x 

n=0 



Theorem 4. For P an arbitrary N -element point set in [0, l] s and w = (w , 
R we have 

(logA0 (s ~ 1)/2 

e(U(K™ )S );P;w) > C(a, 7 ,s)^^ 

where C(a,~y, s) > depends on a,~y and s, but not on N and w. 



wn-i) e 



Proof. We follow the proof of Temlyakov [31, Lemma 3.1]. Let (3 := ^2^=o w n- If /5 = 
then e(%(K™^ s ); P; w) > 1, since for / = 1 the integration error is 1. In this case the 
result holds trivially. Thus we can assume now that /3 ^ 0. 
We have 

e 2 (U(K^ s );P;w) 

. N-l N-l 

= / K(x, y) dx dy - 2 V] w n / K (x, x n ) dx + V] w n w n >K(x n ,x n/ ) 

J V>>V n=0 "V 1 ]* n,n'=0 

^ Wn J]cos(7rA; j a; jin ) 

n=0 j=l 

where a; n = (xi,„, . . . , x s>n ). 

For m = (mi, . . . , m s ) G Nq and |m| := mi + - ■ ■ + m s we will now construct a function 
^(2/) := Sim|=t F m (y)i parametrized by the points x n and weights w n of the arbitrary 
cubature rule, to obtain a lower bound on the worst-case error. For this we will pick its 
cosine coefficients to be bounded above by r a)7iS (fe). Let the integer t be chosen such that 

2N < 2* < AN. 

Let a := [a] + 1 and let / : R — > R be the a-times differentiable function 

V +1 (l -x) a+1 for0<x<l, 







otherwise. 



(8) 
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Note that f(x) > for < x < 1 and supp(/( r >) = (0, 1) for all < r < a. 

For m E N let f m (x) := f(2 m+2 x) and for m = (mi,...,m s ) G Nq and x 
(xi, . . . ,x s ) G R s let 

s 

fm(x) := Y[f mj (Xj). 

3=1 

Then supp(/ m ) = ]lj=i(0, 2" m ^- 2 ). We obtain 



WO) = ^ f /(2^ +2 -) dx = n 2^+2 T /(») = ^ 

where /(/) := J* f(y) dy. For / given by (8) we obtain 

((a+1)!) 2 



WW 



1(f) = B(a + 2,a + 2) 



(2a + 3)! 



where B denotes the beta function. 
For k 7^ we have 



= / f{2 m+2 x)V2cos{irkx)dx 



-m — 2 



-J= J f(2 m+2 x) {e^ kx + e~ nikx ) dx 

f f(y) ( e 2 ^ m ^ + e- 27rifc2_m ^) dy 



2 m+2 v / 2 Jo 
f{k2~ m - 3 ) + f(-k2- m ~ 3 
2 m + 2 V2 



where 



f(h)= I f(x)e- 2 - ihx dx 



o 



denotes the Fourier transform of /. Since, by definition, /M(0) = = for all 

< t < a, and / is a-times differentiable, repeated integration by parts shows that for 
any m G No we have 

|/(A;2- m - 3 )| < C a min(l, {k2~ m ~ 3 )- a ), 
where the constant C a > depends only on a (and /). Thus we have 

\fm,Uk)\ < C a 2^ m - 3/2 min(l,(A;2- m - 3 )- a ) 
< C' a 2- m mm(lX m r a /2,i(k)). 

This bound even holds for / m , C os(0) if C' a is large enough. For the multivariate case we 
have the bound 

s 

|/m,co S (fc)| < C(a,s)l[2-^mm(l,2 a ^r a/2A (k,)) 

i=i 
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s 



= C(a, s) 2 (a - 1)|m| Jj2" am ' min(l, 2 am >r a/2yl (k j )). 

3=1 

By summing |/ m , CO s(fc)| 2 over all choices of m where \m\ = t we obtain 

s 

mSNg mGNg j = l 

ttx|— i m|=t 

s oo 

< 2 2(a - 1) 'C ,2 (a, s) J] S 2_2am ^t 1 ' 22amr a,i(%))- 

j=l m=0 

The last sum can now be estimated by 

oo 

^2- 2Qm min(l,2 2am r ail (A;,)) 

m=0 

2 2{a - a)m r a , 1 (k j ) + ^ 2 ~ 2am 

0<m<(log 2 r-\(kj))/2a m>(log 2 r-j(fc 3 ))/2a 

r^ a , 1 (fc J -)2 2 (— ) - I r a;1 (%)2 2 " 
- 2 2 0~ a ) - 1 a ' 11 3 ^ 2 2a - 1 

2 2a 



< r a>1 fo) (1 + 22Q _ 1 

< 3r a>1 (kj). 

Thus, since 2N < 2* < 47V, we have 

~ 2 2t ~ 

r a , 1 Ak)>C,{a n) s)2- 2 ^ t ^ |/ m , C os(fc)| 2 > d(a, 7 , s)— ^ |/ m ,cos(fc)| 2 . 

-m.| =t |m.|— t 

Now for a; = (xi, . . . , x s ) and y = (y 1 , . . . , y s ) and for u C [s] define 

x(±) u y = (z!,...,z s ), 

where Zj = Xj + yj if j G u and Zj = Xj — yj if j u. Define the functions 

iv-i x 
F m ,u(y) ■= Yl w nfm(x n (±) u y), F m (y) := — ^ F m>u (y), 



2 s 

n=0 uC[s] 



and the sets 



B m ,u-={y^[0,l] s :F m ^y)=0}, 
B m ■= {y e [0, l] s : F TOitt (y) = for all u C [s]} = f| B m>u . 

uC[s] 

Denote with -B^ u the complement with respect to [0, l] s . Then for X s the s-dimensional 
Lebesgue measure we have A s (supp(F mjU )) = A S (S^ IU ). Since supp(/ m (cc n (±) u t/)) as 
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a function of y is contained in the interval JX/eiX - x .?>> — x h n + 2 mj 2 ) Ylje[s]\u( x j^ ~ 
2~ m i~ 2 ,Xj jn ) we have 

AT-1 

snpp(F m , u )C |J J](- x . n ,- x . )n + 2-^- 2 ) J] (x j)n -2-^- 2 ,x j)n ). 
n=o jeu je[s]\u 

Thus A s (supp(F m)U )) = X s (B^ n ) < N2^ m ^ 2s . Now, for all m satisfying |m| = t we 
obtain 



/ 



A s (5 m ) = l-A s (££j = l-A s 



U ] > i - £ iV2Hm| " 2s = 1 " JL > V2, 



2|m| + 

\uC[s] / uC[s] 



since 2N < 2* < 47V. 

We can expand F m (y) — Jj ^ F m (y) dy in terms of the coefficients fm,cos(k): 

1 

^7 £ £ W nfm(x n (±) u f/) - /m,cos(0)/3 
u C[s] n=0 

= E fm,cos(k) £ II ^(^J 'fen + %)) II <Xx(Trkj(x j>n - Vj)) 

n=o fceNg\{o} uc[ s ]jeu je[s]\u 

£ /mcos(fc)(v^) |fc| ° f[ COS ( 7rfc j( a: J> ~ y Jl + COS (^fe» + Vj)) 
n=0 fceN«\{0} j=l 
N-l s 

= E w ™ £ /m,cos(fc)(v / 2) |fc| ° JJ(cos(7rA; i a; iin )cos(7rA; j i/ j )) 

n=0 fc£N»\{0} j=l 

(AT-1 s \ s 

^w^cos^-i^) J JJcos(7rA; i y i ). 

n=0 j=l / i=l 

Thus, by definition of B m , we have 

A s (B m )|/ m ,cos(0)| 2 /3 2 = J b ^EE^^W^)-/'"^ d ^ 



uC[s] n=0 

- / W nfm{x n (±) u y) - 7 m , cos (0)/3 

-V 1 ]* V «<=[«] n=0 

(AT-1 s 
£ »n JJcos(vrA: 
n=0 j=l 



dy 



v j x j,n) 



We are now ready to piece this all together to obtain 

2 



fceNg\{0} \n=0 j=l 
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N-1 s x " 



2 2t ~ 



N 2c 

"iSNg fceNg\{0} \n=0 j=l 

\m\=t 



> 



2 2t ~ 
[1-P) 2 + Cl (a,j,s)— K(B m )\f m , cos (0)\ 2 P 2 







| m \=i 

> (1 - /3) 2 + C7 2 (a, 7, s)f (I(f)) 2s j^l~ 2t -' s £ 1 

> (l-/3) 2 + C 3 (a,7,s)/3 2 iV- 2a/f 



|7Tl|=t 



t + s - 1 
s - 1 

Set A := C 3 (a, 7, s)A r_2a (*^ 1 ). Then the last expression can be written as (l-(3) 2 +A(3 2 , 
which satisfies 

e\nK™ s )- P- w) > (1 - /3) 2 + A/3 2 > > C 4 (a, 7, S )iV" 2 " ^ + 8 _~ 1 

which implies the result, since t > log 2 (A^). □ 

5 Numerical results 

In this section we show some numerical examples of applying lattice rules, tent-transformed 
lattice rules and symmetrized lattice rules to some test functions. For this we use the 
lattice sequence from [16] which was constructed to give 3rd order convergence in a Ko- 
robov space. It is a 10-dimensional base 2 sequence with a maximum of 2 20 points and is 
comprised of embedded lattice rules with sizes 2 m for m = 0, . . . , 20. This lattice sequence 
was also used in [27] for some higher order convergence tests. 
We report on two test functions to show some effects: 

g S;W (x) := f[ (l + ^ (-10 + 42x 2 - 42x 5 j + 21x 6 ^ , 

KA X ) '■= II ( l + Y ( 31 - Mx ) + 8x i + 70x j - 2&r j + 8x ] - 16 - 16 sin(x i )) N ) . 

Both functions integrate to 1 over [0, l] s . The parameter w acts like a product weight wK 
All tests use 2 20 (plus 1 for the symmetrized rule) function evaluations for their final 
result. In Figure 1 we report the actual running time in microseconds of optimized C++ 
code (to accommodate for the difference in which integration nodes are constructed). 
Every mark represents an approximation with 2 m function evaluations or 2 s_1 2 m + 1 
for the symmetrized rule. From this it can be seen that all three methods approximately 
have the same cost in this implementation, but in general this is dependent on the relative 
differences in the time for generating lattice sequence points, symmetrization and function 
evaluation. For completeness we note that the tests were run on a 1.8 GHz Intel Core i7 
and compiled with the clang++ 3.1 C++ compiler. Tests on an older 2 GHz Intel Xeon 
compiled with the g++ 4.8 compiler gave similar looking results (but slower). The C++ 
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Figure 1: Left column: function g SjW (x) in 8 dimensions for w = 0.1, w = 0.9 and w = 2. 
Right column: function h StW (x) in 10 dimensions for w — 0.1, w — 0.9 and w — 1. 



source code and the raw data (for more test functions than shown here) can be downloaded 
from the KU Leuven Lirias document repository. 

The function g S)W {x) is a product of 1 (B 6 (xj) + E 5 (xj)), where B 6 is the degree 6 
Bernoulli polynomial and E 5 is the degree 5 Euler polynomial. From their Fourier ex- 
pansions [1] follows that we expect 3rd order convergence for both the tent-transformed 
and symmetrized lattice sequence, while we only expect 1st order convergence for the 
standard lattice sequence. We show some results in the left hand column of Figure 1 for 
s = 8. When the weight parameter is close to 1 (w = 0.9 in the middle panel) we notice 
that the performance of the symmetrized and the tent-transformed rules is similar. If w 
is much smaller (w = 0.1 in the top panel) then the tent-transformed rule wins, while if 
w is larger (w = 2 in the bottom panel) then the symmetrized rule wins; these effects are 
more pronounced when the dimension gets larger. 

The function h StW (x) is a more arbitrary function (including sin(xj) and E 7 (xj)) for 
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which its cosine expansion converges only like k~ 2 . We thus expect 1st order convergence 
for the tent-transformed and symmetrized methods. For s — 10 the behavior is illustrated 
in the right hand column of Figure 1. Here again it can be seen that small values of w 
give the advantage to the tent-transformed lattice rule, see the top panel where w = 0.1. 
However, already for w = 0.9 the symmetrized rule outperforms both other methods. 
The bottom panel shows the case w = 1 where it can be seen that the number of points 
needed to reduce the error below 1 is already getting quite high. In such a case there 
is an exponential dependence on the number of dimensions and the symmetrized rule is 
then the preferred choice if the dimension is small enough (this can also be seen in the 
left bottom panel which has w = 2 for g s , w ( x ) in s = 8 dimensions). We remark that 
the graphs, where one replaces the time in the abscissa by the number of points N in 
Figure 1, are very similar to the ones shown in Figure 1. 

6 Conclusion 

In this paper we have shown that it is possible to obtain a convergence rate of N~ a+S 
for any 5 > for sufficiently smooth integrands for lattice-type rules also for nonperiodic 
functions. Previously, the only QMC rules with such properties where higher order digital 
nets [9] . However, the function spaces we consider here are smaller than the usual smooth 
Sobolev spaces used for higher order digital nets. 

Since for smoothness 1 the unanchored Sobolev space and the half-period cosine space 
coincide, we obtain that tent-transformed lattice rules achieve the same convergence be- 
havior and tractability results as lattice rules in Korobov spaces. In contrast to previous 
results this technique does not need randomization. 
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